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We studied the diluted magnetic semiconductor by the self-consistent Green's fimction approach, 
which treats the spin-wave kinematics appropriately at finite temperatures. Our approach leads to 
a simple formula for the critical temperature in a wide range of parameter space. In addition, the 
magnetization curve versus temperature in some regimes is concave, which is dramatically different 
from the usual convex shape. Finally, we discuss the possibility of generalizing the current theory 
to include the realistic band structure, electronic correlations and disorders in a systematic way. 

PACS numbers: 75.30.Ds, 75.40.Gb, 75.50.Dd 



Diluted magnetic semiconductor (DMS) has attracted 
intense attentionsj^ 21 for its potential applications in 
spintronics devices. Robust ferromagnetic order has been 
observed in (Gai_a;Mn;j;)As up to 110 KjHi]. Magnet- 
ically doped wide bandgap semiconductors and oxides 
such as GaN, ZnOjTi02 even exhibit ferromagnetism at 
room temperature|5|,|y,t3, although the magnetization is 
less robust than the doped III-V semiconductors. While 
lots of efforts are focused on the search of optimal ma- 
terials virith enhanced critical temperatures Isl, it remains 
a challenging task to describe the coupled localized mo- 
ments and the itinerant carriers in an appropriate way. 



Part of the difficulty lies in the fact that the spatial 
fluctuations of the ferromagnetic order is large at finite 
temperature and the usual Weiss mean-field theory does 
not work. To account for th e sp atial fluctuations, nu- 
merical approaches^, 0, 0|, such as local density 
functional approximation and Monte Carlo simulations, 
are quite helpful in estimation of various thermodynam- 
ical properties. However, it is rather difficult to study 
the electronic transport by these numerical approaches. 
Analytical approaches[l^ [ij, provide partial de- 

scriptions in several particular limits by treating the im- 
purity spin semiclassically or ignoring the interactions 
and the kinematic constraints among spin waves. At fi- 
nite temperatures, the average spin- wave density is large 
so that these approximations are no longer appropriate. 
It is therefore desirable to develop a spin-wave theory 
which works at finite temperature. 

In this Letter, we adapt the self-consistent Green's 
function approach to describe the fluctuating spin cor- 
relations at finite temperature. Since the kinematic con- 
straint of spin waves are treated exactly in the equa- 
tion of motion, this method can be applied to a wide 
range of parameter space, even very close to the critical 
temperature^, 0|. Indeed, the critical temperature is 



determined by the simple formula, 
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where c is the density of impurity spins. The spin- wave 
dispersion Qc{p) and the impurity spin polarization etc 
are both determined self-consistently near the critical 
temperature T ^ Tc. The trends of the critical tem- 
perature upon the change of parameters are studied in 
detail later. 

We model the DMS by the Hamiltonian, containing 
only the kinetic energy of itinerant carriers and the ex- 
change interaction between the itinerant and the local- 
ized impurity spins, 



H = Ho + J / d^r S{r) ■ s{r), 



(2) 



where J > is the strength of exchange interaction. The 
impurity spin density is S{r) — (5^^-'(r — Ri)Si while 
the itinerant spin density is s(r) = ^^p^r){a■ /2)il!{r). The 
band structure of the itinerant carriers is described by 
Hq, which depends on the host semiconductors. Since 
our emphasis here is how to cope with spatially fluctu- 
ations appropriately, the dispersion is taken as the sim- 
plest parabolic band, Hq — /2m* . Generalization to 
more realistic but complex band structures, such as the 
Luttinger model, can be achieved straightforwardly. 

Since the impurity spins are randomly doped into the 
host semiconductor, their positions are random. If the 
disorder is strong, the itinerant electrons are localized 
and the percolation approach 0| would be more ap- 
propriate. However, we are interested in the metallic 
regime where itinerant carriers are delocalized. The dis- 
order also plays a crucial role in smoothing out the im- 
purity spin density S{r). For instance, the magnitude of 
the spin density is smeared after coarse-graining. 
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FIG. 1: The magnon excitation spectrum at zero tempera- 
ture. The impurity density c — 1.0/nm^ and c* — 0.01/nm''. 
The shaded area represents the regime where the imaginary 
part of the self energy is not zero but negligibly small. 



where c is the impurity spin density and the average 
is taken over random locations of the impurity spins. 
Therefore, the presence of weak disorder allows a field- 
theory description in the continuous limit. 

The dynamics of the impurity spins is described by the 
thermal Green's function, 



Dir,T) = ((5+(r,T);5-(0,0))) 

^ -e(r)([5+(r,T),^-(0,0)]). 



(4) 



The equation of motion for -D(r, r) would involve more 
Green's functions of higher orders. The exact solution 
then involves the Green's functions of all orders and 
is not feasible in general. However, within mean-field 
approximation, the higher-order Green's functions can 
be decomposed into simpler ones and eventually a self- 
consistent solution is possible. For the spin-wave propa- 
gator D(q,ii>n) in momentum space, the mean-field de- 
composition simplifies the equation of motion, 

J{Sz) 

c* V ^ — ' 

k 

where c* is the density of itinerant carriers. Notice that 
it only involves one additional Green's function F{k, k + 
q,ivn) = {{'ip\ik)ipi{k+q); S~{0,0))). Applying the same 
trick again, the equation of motion for F{k, k + q, ivn) is 



F{k, k + q, ivn) 



Jc* /T(fc)-/i(fc+P) 
2 zj/„ -f A -I- efc - Efc+p 



D{q, ij/„). 



(6) 

where f^^i{k) = [eP(^k±A/2-f,) _^ ^-^-i jg ^j^g Feimi dis- 
tribution for itinerant carriers with different spins. The 
Zeeman gap A = J{Sz) in the electronic band structure 
is due to the ferromagnetic order of the impurity spin 
and has to be determined self-consistently. 



From Eqs. (0) and ©, we can solve for the spin- wave 
propagator. Upon the Wick rotation, iun ft + irj, the 
spin-wave dispersion is identified as the simple pole in 
the Greens' function, 



n - j{sz 



2V ^n + A + ek- £k+p + iv 



J A ^ 



0. (7) 



The average spin densities (s^) and (S^) in the above 
equation remain unknown and need to be determined 
self-consistently. 

The average itinerant spin density (s^) is just the dif- 
ference of the spin densities in majority and minority 
bands, split by the Zeeman gap A. The relation between 
the average localized spin density (Sz) and the spin- wave 
dispersion fl(k) is more subtle due to the non-trivial spin 
kinematic constraint. From Callen's formula. |l9j 

l/c,v_c / (25-|-l)K^)^^+i 
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Here {n^^) = (1/cV) X^fcl^^"^*^ - is the average 
number of spin waves over all momenta. At low tempera- 
tures, the average spin- wave number is small and the last 
term in Eq. ((SJ can be safely ignored. The magnetization 
show the T^/"^ behavior as in the independent spin-wave 
approximation. As the temperature approaches the crit- 
ical regime, the spin-wave density becomes large and the 
kinematic constraint becomes important. Solving Eqs. 
{Tj) and ©J both the spin- wave dispersion n(k) and the 
averaged impurity spin density {Sz) are obtained self- 
consistently. 

The dispersion from Eq. ^ has two branches due to 
the presence of both itinerant and localized spins. Be- 
cause the gapless spin-wave fluctuations dominates, the 
optical branch can be safely ignored. Besides, the spec- 
tral weight of the optical mode is also small due to the 
dilute density of itinerant carriers. In Fig. the lower 
branch of spin-wave dispersion at zero temperature is 
shown, corresponding to the Goldstone excitation of the 
ferromagnetic order. Our numerical results show that 
the imaginary part of the spin-wave self energy, due to 
the presence of Stoner continuum (the optical branch), 
is negligible in all regimes and the dispersion is sharp. 
Thus, it is a reasonable approximation to assume the 
spectral weight is carried by the gapless spin-wave exci- 
tations only. 

The magnetization curves at different spin densities 
c*, c are plotted in Figs. |21and|21 Throughout this Letter, 
the exchange coupling and the effective mass are fixed at 
typical values J = 0.15 eV nm^ and m* — 0.5mey . First 
of all, we study the temperature dependence of the mag- 
netization at different impurity spin densities c, while 
the ratio of itinerant and localized spin densities fixed 
at c*/c = 0.1. It is quite interesting to notice that the 
magnetization drops more dramatically near the critical 
regime when the impurity spin density is large. As shown 
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FIG. 2: Temperature dependence of spin polarizations at 
different impurity densities c. The ratio of itinerant and im- 
purity spin densities are fixed at c*/c = 0.1 for all curves. 
The inserted figure shows the monotonically increasing criti- 
cal temperature Tc at different densities. 
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FIG. 3: Magnetization curves at different itinerant spin den- 
sities c* , while the impurity spin density is fixed at c = 1 
nm~^. The critical temperature Tc, as shown in the inset, 
reaches the maximum at an optimal concentration of the itin- 
erant spin density. 

in Fig. [3 the critical temperature, computed by Eq.^ in- 
creases monotonically as the densities increase ,20] . This 
monotonic increase in the critical temperature is qualita- 
tively (not quantitatively) the same as the Weiss mean- 
field theory and also agree with the experiments. 

On the other hand, if the impurity spin density is hold 
constant, say c—\ nm~^, both the magnetization curve 
and the trend of the Curie temperatures show interest- 
ing behaviors, which deviate from the Weiss mean-field 
theory. Starting from extremely dilute density of the 
itinerant carriers, the magnetization curve is concave as 
shown in Fig. |3 in contrast to the usual convex cur- 
vature in Weiss mean-field theory. This concave shape 
resembles the magnetization curve in percolation theory 



near the boundary of the metal-insulator transition and 
is often used to be an indication of localization in the 
presence of disorder [H However, our results show 
that a concave magnetization curve is not necessarily tied 
up to the localization tendency. We emphasize that the 
disorder of impurity spin is included in our approach only 
through the coarse-graining procedure and the transport 
of the itinerant carriers is assumed ballistic here. So it is 
rather surprising, by including the spin- wave fluctuations 
appropriately at finite temperatures, the magnetization 
curve is concave at dilute densities. 

In the regime where the magnetization curve is con- 
cave, the critical temperature increases with the carrier 
concentration c*. Beyond an optimal density, c* ~ 0.01 
nm^"^ (the impurity spin density is fixed at c = 1 nm^'^), 
the critical temperature reaches the maximum and starts 
to fall back as in the inset of Fig. 13 The suppression of 
the critical temperature is mainly due to the oscillatory 
RKKY interaction at higher concentrations. The effec- 
tive coupling between impurity spins is frustrated and 
not longer purely ferromagnetic. The existence of an op- 
timal density is totally missed in the Weiss mean-field 
theory, where both the spatially varying fluctuations and 
the quantum frustrations are ignored. 

Beyond the optimal density, not only the critical tem- 
perature falls back, the shape of the magnetization also 
undergoes an interesting change. There are two signifi- 
cant features. One is the magnetization curve becomes 
convex again. Another feature is that the magnetization 
has a very sharp decrease near the critical temperature, 
making it look almost like the first-order phase transi- 
tion. To make sure that the sharp drop is not an arti- 
fact of numerical errors, the data points near this regime 
are chosen very closely to ensure we still have enough 
points in the steep regime. In particular, for c* = 0.05, it 
is spectacular that the magnetization drops 50% within 
IK. It is not clear at this point what is the underlying 
mechanism behind this dramatic suppression of magne- 
tization. A more sophisticated theory starting from the 
critical point might be able to address this interesting 
dive of magnetization near the critical temperature. 

Since the global trend of the critical temperature is 
important, it is plotted in Fig.^ Comparing with Fig. 2 
in Ref.j^, the self-consistent Green's function approach 
produces a more complicated landscape. In the Weiss 
mean-filed theory, the critical temperature is 

^MF ^ XP S{S+l)J^N 

(5>b/2)2 12 ' ^ ^ 

where g* is the (7-factor of the carriers and xp is their 
Pauli susceptibility, which is proportional to the effec- 
tive band mass. The above formula would produce a 
profile with monotonically increasing Tc as the densities 
become larger. It is clear from Fig. 21 that the profile 
of the Curie temperature has a ridge, roughly along the 
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FIG. 4: The trend of the critical temperatures at different 
densities. 

curve c ~ 100c* with the given parameters in this Let- 
ter, and is qualitatively different from previous studies. A 
closer check would find that our approach also produces 
quantitative differences. Comparing with the estimates 
of Curie temperatures, previously studied by one of the 
authors in Refs. ^3 and 0|, the present approach gives 
a lower Tc in all densities regimes. So it seems that the 
inclusion of spin kinematics at finite temperatures is not 
only important to get the right trend of the critical tem- 
peratures, but also crucial in estimating their values. 

Finally, we address the important aspects of physics 
which are left out so far. To make quantitative com- 
parison with the experiments, it is crucial to adapt the 
realistic band structure of the host semiconductors, for 
instance, the six-band Luttinger model for GaAs. The 
inclusion of the more complex band structure would not 
cause formidable messes in calculating the Green's func- 
tion self-consistently. While most of the conclusions 
drawn from the simple parabolic band should remain 
valid, a more realistic band structure is desirable for mak- 
ing quantitative predictions. The electronic correlation 
can also be included at the mean-field level in a system- 
atic way. Since Coulomb repulsive interaction stabilizes 
the ferromagnetic phase, the Curie temperature is ex- 
pected to be higher. To include the disorder is more 
subtle. In the diffusive regime, one can replace the bal- 
listic propagator of the itinerant carriers by the diffusive 
one to account for the disorder effects. 

In conclusion, we employ the self-consistent Green's 
function approach to study the DMS and derive a gen- 



eral formula for the Curie temperature. In addition, we 
demonstrate the interesting crossover of the magnetiza- 
tion curve from concave to convex, which is not driven 
by the localization effects. 
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